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The amplitude and phase are calculated for oblique propagation across a coastline with 
a sloping beach. In this case, the land and sea are taken to be plane surfaces and the beach 
slope is constant. It is shown that the reflected wave may be quite significant and it 
has a fundamentally different character from the reflected wave in the case of a flat-lying 
coastline. 



1. Introduction 

In part I [Wait, 1963] of this series, the propagation 
of radiowaves across a flat-lying coastline was con- 
sidered in some detail. It is the purpose of this 
paper to extend the analysis to include the effect of a 
gradual elevation change between land and sea. 
Numerical results are given for a beach which has a 
plane or constant slope. The same problem has 
been treated by Feinberg [1946] and more recently 
by Kalinin [1958]. Our results, derived in a different 
manner, appear to be in agreement with the Soviet 
work. 

2. Formulation 

The general situation is illustrated in figure 1. 
With respect to a Cartesian coordinate system, the 
xy plane is taken to be the plane surface of the sea to 
the right of the waterline at x=d . The elevation of 
the land is then defined by 



and 



z=z for £<0, 
z=z (d —x)ldo for Q<^x<^d . 



To simplify the discussion, it is assumed that the 
electrical characteristics of the land may be described 
in terms of surface impedance Z right up to the water- 
line atx=d - The medium to the right (i.e., x>rf ), 
which is the sea, is described by a surface impedance 
Z'. As in part I, the transmitter at A with coordi- 
nates (—x , yo) is regarded as a vertical electric dipole, 
of effective height h a , on the surface of the land. 
The receiving antenna, of effective height h b) is 
located at B with coordinates (d i} 0) where d x may 
be positive or negative. 

It is apparent that a rigorous solution of the 
problem would be extremely difficult. Furthermore, 
even if it were available, the idealization of the model 
would limit its usefulness. Therefore, an approxi- 
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mate approach is adopted which leads to a relatively 
simple formula for computation. 

Provided the slope of the beach is small, it can be 
expected, on the basis of physical consideration, that 
the main influence of the elevation change is to 
modify the tilt of the electric field. For example, if 
the slope in the x direction is defined by y x , then the 
modification of the horizontal electric field E x is 
approximately equal to —y x E z . To within a first 
order, the horizontal electric field E y is not changed 
if the slope y y is zero. 

In order to achieve a further simplification, the 
source dipole is assumed to be sufficiently removed 
(to the left in fig. 1) that the incident wave is nearly 
plane in the vicinity of the coastline. Thus, locally 
the incident field is proportional to exp [ikSiy—ikCix] 
where C x and Si are the cosine and sine of the angle 



A(-x ,y ) 



> 



B(di,o) 



CLAND) 
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Figure 1. Plan and side view of the mixed path showing the 
location of diyoles A and B. 



The scales are exaggerated 
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of incidence 6 at the coastline. It is then a simple 
matter to describe the tilting of the electric field in 
terms of the modification of the effective surface 
impedance. Thus, over the plane z=z and in the 
interval 0<jr<Aj it is found that the impedance 
boundary conditions may be written 



where 
and 



-Z x H y and E y =Z y H x , 



Z x =Z—y x rio/Ci, 



Zy Z. 



Consequently, we have, in effect, replaced the 
sloping beach by a horizontal plane surface with a 
modified surface impedance. Boundary conditions 
of this type can be deduced directly from the results 
of Feinberg [1946] after a long and involved deriva- 
tion. According to Feinberg, the restrictions on 
the slope are equivalent to requiring that y x <^<C\ 
and z <^<^^d \/27r. 

3. Integral Representation 

The mutual impedance between the dipoles A 
and B is again denoted by z m if the whole ground 
plane were flat and had a constant surface impedance 
Z. The change of the mutual impedance which 
results in the differing electrical characteristics of 
the sea is denoted Az m . In part I, the calculation 
of Az m was carried out under the assumption of a 
sharp boundary or sudden change from Z to Z / at 
x=0. The change of the mutual impedance result- 
ing from the elevation change is denoted 8z m . In 
this case, 8z m is expressed as a surface integral over 
the strip 0^x^d . Thus eq (1) of part I is now 
replaced by 



te. 



•4JJV 

strip 



-yzVo/CJHavHiydx dy, 



(2) 



where H ay is the y component of the tangential 
magnetic field over the strip if the surface were 
unperturbed and H' bv is the tangential magnetic 
field over the strip under perturbed conditions. In 
formulating this integral, the current at the terminals 
of dipoles A and B is /. 

The simplification and reduction of the above 
representation for 8z m is carried out in the manner 
described in part I. There is an essential difference 
here in that only the y components of the magnetic 
field are involved. Thus, omitting numerous details, 
it is found that 



8z 

z 



XH^ikl^-xlCJdx. (3) 



Actually, this result is valid even when y x is some 
smooth function of x. However, it is necessary that 
y x have a magnitude that is small compared with 
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Figure 2. Real and imaginary parts of the function F(ai) as 
a function of ai for clq = 0.6~. 

unity corresponding to a very gentle slope. For 
the problem described here, 



lx=—j [u(x)—u(x—d )l 



and 



<>y x 
bx 



2o 



[8(x)-5(x-d )], 



(4a) 



(4b) 



where u{x) and 8(x) are the unit step and unit impulse 
functions, respectively. The integrations, with re- 
spect to x, may now be carried out readily if the 
identity, 

~ [ae**°(H?>(a) TiH^(a))] = e^m 2) (a), (5) 



introduced in part I, is again utilized. Thus 

-(i) FW ' 



Q~Z m _ 
Zm 



where 

F(a 1 )=J(a 1 )—f(a 1 —a ) for a_>a , 

=jf(a_) + (j{a Q —a l ) for 0<a!<^ , 

=g(ao— «_)— g(— ai) for«!<0, 
where a =&CWo an( i oti = JcCidi. 



(6) 

(7a) 
(7b) 
(7c) 
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Figure 3. Real and imaginary 'parts of the function F(ai) as 
a function of on for clq=1.0. 

The basic functions/ and g arc defined by 

jf(x) =xe ix [m 2) (x) -illi 2) (x)]+ie ix H!> 2) (x) , (8a) 
and 

(fix) =xe~ ix [m 2) (x) + ill[ 2) (x))-ie- ix m 2) (x) , (8b) 

where % is the general argument (which is positive 
real). 2 

The function F(ai) is proportional to the fractional 
change of the field at B for a source dipole at A. 
Conversely, F{a x ) may be regarded as the fractional 
change of the field at A due to a source at B. The 
reciprocity is an inherent feature of a mutual im- 
pedance formulation. It is emphasized that the 
results are valid only when \8Zm/'z m \«l. 

The parameter a Q is proportional to the width of 
the coastal strip whereas a y is proportional to the 
distance d x which is measured from x=0. The real 
and imaginary parts of the function F{a x ) are shown 
plotted in figures 2 to 6 for a range of values of a 
from 0.5 to 5.0. The curves all have an oscillatory 



2 It may be noted that for x<<1, 

f( x )^-g(x)^Yo(x)-xYi(x)^-- [log. (2/ x ) -1.5772]. 
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Figure 4. Real and imaginary parts of the function F(ai) as 
a function of ai for a u = 2.0. 

behavior for negative values of a x . This is a mani- 
festation of the interference between the incident 
wave and the wave reflected from the coastline. 
For positive values of a u corresponding to points 
beyond the coastline, F{a x ) has a monotonic 
character. 

As may be observed in figures 2a, 3a, 4a, 5a, and 
6a, the real parts of the function F(a x ) exhibit 
singularities at the edges of the coastal strip (i.e., 
at «! = and ai = a ). Actually, the results are not 
valid in the vicinity of these singular points since 
\oz m /z m \ is then no longer small compared with unity. 
Nevertheless, it can be expected that the real part 
of the field function F{a { ) should show marked 
changes in these regions. 

The phase of the function F(ai), as indicated in 
figures 2b, 3b, 4b, 5b, and 6b, is a continuous func- 
tion right across the coastal strip. It is particu- 
larly interesting to note that the pronounced ampli- 
tude changes at the edges of the coastal strip are not 
accompanied by marked phase variations. Of 
course, this fact is compatible with the behavior of 
the rigorous field solution in the vicinity of the apex 
of a perfectly conducting wedge [Wait, 1959]. 
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Figure 5. Real and imaginary parts of the function F(ai) as 
a function of ai for a = 3.0. 

It is of practical importance to observe that the 
phase perturbation resulting from the coastline is 
reduced to negligible proportions within several 
wavelengths of the coastline. This is in sharp dis- 
tinction to the phase variations which result from 
the conductivity contrast between land and sea. 
As indicated in part I, the phase beyond the coast- 
line (away from the transmitter) increases contin- 
uously with distance. 

In order to summarize conveniently the field 
variation behind the coastline, the real and imaginary 
parts of F(a ] ) 1 for positive values of a u are plotted 
in figures 7 a and 7b. 

To facilitate further calculations for the influence 
of the sloping beach, the general functions /(%) and 
g(x) are presented in graphical form in figures 8 a 
and 8b. These functions, defined by eqs (8a) and 
(8b), are now written in the form 



y(x)=/,(x)+«y«(x), 



(9) 



and 



g(x)=gr(x)+igi(x), (10) 

where f r , f t , g T} and g t are real functions of x- As 
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Figure 6. Real and imaginary parts of the function F(ai) as 
a function of ai for ao= 5.0. 

indicated, the function g(x) is generally oscillatory. 

4. Refraction Effects Resulting From the 
Sloping Beach 

It is interesting to note that the sloping beach 
actually causes a refraction error which is more or less 
additive to the error resulting from the conductivity 
contrast with the sea. As in part I, it is assumed 
that the field incident on the boundary at x=0, has 
the form 

E =e- ikClX e ikSiV , 

which is appropriate if the transmitter at A is suffi- 
ciently far removed to the left in figure 1. Under 
the further assumption that the total perturbation 
of the field E is small, it follows readily that 
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Figure 7. Real and imaginary parts of F(ai) for positive 
values of (ai — a ) for various values of a . 



where Az m /z m is the relative perturbation resulting 
from the conductivity contrast and 8z m /z m is the 
relative perturbation resulting from the sloping beach. 
Following the line of reasoning given in part I, the 
total refraction or bearing error 86 T may be expressed 
in the form 

86 T =8S+86', (12) 

where 86 is given explicitly by eq (36) of part I, while 
WeffiCisr- bn(—\ (13) 

Oai \ z m J 

Restricting attention to the region beyond the coast- 
line (away from the transmitter) and employing eq 
(6) it follows, without difficulty, that 



86'J^± (|) [sin a 1 Y 1 (a 1 )+COS aJM) 

—sin (ai--ab)Fi(ai— <xq) — COS (a x — oq) Ji{a x — oj )], 

(14) 

where the symbols have their usual meaning. In 
particular, J x and Y x are the Bessel and Neumann 
functions, respectively, of order one. Provided 
{a\ — a )>^>l) the above result may be simplified to 

W J^(*J\\ i LI 

: 2 WLk(ai"«b)]* [«ft]*J 

sin o (cos q q )\ / Zq 



2tt* 



wL[*(di-do)]* [Mi] J* (15) 



If, in addition, di^>y>d corresponding to observa- 
tions sufficiently far from the coastline, the above 
formula simplifies, even further, to 



86'^ 



sin 6 (cos O )* z 1 



47P 



d\ (kdi) 1 * 



(16) 



which is independent of c/ . 



It is interesting to note that 86' varies approxi- 
mately as {di)~' ir \ whereas 86 varies approximately 
as (di)~ 1/2 . Thus, it is concluded that the bearing 
error resulting from topographical features at the 
coastline is probably of minor consequence in com- 
parison with the influence of conductivity contrast. 
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